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I. INTRODUCTION 

The gravitational theory with higher derivative terms [1, 2] is one of the most interest¬ 
ing subjects to study in physics for many years. The treatment of Einstein gravity with 
perturbative held quantization naturally leads to inclusion of such nonlinear terms in the 
spacetime curvatures in effective actions. 

The generalization of Einstein gravity has been developed extensively by many authors. 
The Lanczos-Lovelock gravity [3-6]^, which is dictated by the action with the dimensionally 
continued Euler forms, is dehned in a generic dimensional spacetime and it describes a 
massless spin-two huctuation despite of the existence of higher order terms in its action. In 
a few decades ago, it is discovered that the effective held theory arising from string theory 
involves similar terms at least the lowest-order correction to the Einstein-Hilbert term [8, 9]. 

The other specihc extension is found in the study on so-called critical gravity [10-12] 
(although the word “critical” originally implied the emergence of Log gravity [13-15]). In 
such a theory of gravity, the behavior of spin-two modes is described by the Lagrangian of 
special combination of higher order terms in the Ricci tensor and the scalar curvature. There 
is a ghost held but no spin-zero mode in the theory. In the three dimensional spacetime, the 
model governed by a similar Lagrangian expresses the new massive gravity [16, 17]. 

In this paper, we consider incorporation of Maxwell helds to the (next-to) critical grav¬ 
ity. The critical gravity has been studied under the assumption of the maximally symmetric 
background spacetime. The higher derivative theory on the other type of background helds 
is worth studying. Inclusion of a non-vanishing hux held brings about the partial compacti- 
hcation with a maximally symmetric extra space and might provide an interesting scenario 
for a very early stage of the universe. 

The assumption of the partially symmetric solution can yield another possibility. It is 
feasible to attach the other type of higher order terms to the Lagrangian which realizes 
the spontaneous compactihcation. Because of the symmetry, use of the Riemann tensor 
is possible in the additional terms. They are expressed by special combinations of Gauss- 
Bonnet terms and Horndeski’s generalized Maxwell term [18]. This can induce a second 
order held equation displaying the cosmological evolution of the background geometry.^ 

^ For a review, see [7] 

^ The applications of the Horndeski’s vector-tensor theory have been discussed in four dimensional [19, 20] 
and also in six dimensional cosmology [21]. 
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This paper is organized as follows. In Sec. II, we build the Lagrangian of our hrst 
model, which involves the quadratic terms of the Ricci tensor, the Ricci scalar and the two- 
form flux held, in a similar way revealed in the study of critical gravity. The spontaneous 
compactihcation in the hrst model is studied in Sec. Ill and cosmological solutions are 
examined in Sec. IV. In Sec. V, we present the second model containing the quadratic term 
of the Riemann tensor. We also investigate the cosmological solution qualitatively in this 
model. The last section is devoted to summary and prospects. 

Throughout the present paper, we deal with our models at the classical level. 


II. A CRITICAL MODIFICATION OF EINSTEIN-MAXWELL THEORY 


In the work on critical gravity, it is found that the Lagrangian of the specihc combination 
of curvatures yields a spin-two massless graviton and a ghost mode; there appears no scalar 
mode as in Einstein gravity. In this section, we propose a model of higher order extension 
of Maxwell-Einstein theory. We will follow the method for constructing the Lagrangian 
examined in our previous work [22, 23]. 

First, we start with the action for the Einstein gravity and electromagnetism in D di¬ 
mensions: 

/o = / d^Xy/^L, = jd^x^[R- FmnF^^ - a] , (2.1) 

where R and A denote the scalar curvature obtained from the metric qmn and the cosmo¬ 
logical constant, respectively. The held strength Fmn is dehned by Fmn = Qm^n — Qn^m 
with a vector held Am and = FmnF^^. The suffices M,N,... vary over 0,1,..., F — 1. 

Applying the variational principle to the action results in the following equations of 
motion: 

Tmn = 0 (2.2) 



(2.3) 

(2.4) 


Next, we construct a Lagrangian of our model as follows: 

L = aL(,-/3r""A„„pQr'’«, 


(2.5) 
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where 

^MNPQ = -^{gMpgNQ + gMQgNp) — _ ^ gMNgpQ ■ ( 2 . 6 ) 

Note that the solution of Tmn = 0 is the solution for the equation of motion derived from 
L, provided that the held equation for the vector held is satished (though the stability is 
not necessarily guaranteed).^ 

The Lagrangian L can then be written as 

L = a[R- - A 

p D jpMN ^ d2 a 1^2 jpMN I D 2, 2 

- P UMN-tt - 4('£) _ 1)'^ ~ + 2(£) _ 1)'^ ^ 

, ^ Tp2 Tp2MN T* + 8 2^2 

Incidentally, we can pick up some ‘critical’ cases: If a = the cosmological 

constant is zero in the D dimensional spacetime; If o; = Einstein-Hilbert 

term (i?) is absent from the D dimensional action; If a = pure Maxwell term 

(F^) is absent from the D dimensional action. 

We should note that the auxiliary held method [14, 16, 17, 23], which can be used in 
critical higher order gravity, also leads to the additional terms, such as 

L' = aLo - /3{2 TmnS^^ - SmnS^^ + S ^), (2.8) 

where is an auxiliary symmetric tensor held and S = S^^qmn- Varying the auxiliary 
held , we obtain 

Tmn = Smn — Sqmn , (2.9) 

or, solving that with respect to Smn, we get 

Smn = Tmn — —-TgMN = ^mnpqT^^ . (2.10) 

Therefore, one can see the equivalence between the Lagrangians L and L'. We can derive 
the linearized held equation using the action I' = j d^xL'. 


^ Generally speaking, we can add the term oc JmJ^ to our Lagrangian, but we restrict ourselves to the 
case without such a term in the present paper. 
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Now, we investigate the behavior of linear fluctuations around the background geometry, 
which satishes = 0. We decompose the metric as follows: 

Qmn = Qmn + hMN ■ ( 2 . 11 ) 

The indices are raised and lowered by the background metric Qmn- Then the trace of the 
fluctuation is expressed as 

h = g^^hMN- ( 2 - 12 ) 

For conveniences, we write down the curvature tensors up to the linear order in Hmn here 

|24]: 

dMTV _ ryMN ^ ^ML iN ^ ^NL iM 

PQ — ^ PQ - 2^ PQi^L ~ 2^ PQi^L 

+l(_VpV^hg + , (2.13) 

jdM _ r)M 3MP r^Q , ^ f}M vL ^ f)L r^M 

— Kn - n NQPp + 2 ^l ~ 2^n^l 

+ , (2.14) 

R = R - R^^hMN + VMVNh^^ -V^h, (2.15) 

where the nabla stands for the covariant derivative and the barred symbols show that they 
are constructed from the background metric Qmn- 

Assuming the background helds satishes Tmn = Rmn — \RgMN — 2 {fIin ~ \F‘^9mn) + 
^AguN = 0, we can use 

5{2S^^Tpq - SpqS^^ + ^2) = 2S^^6Tpq 

+ (higher order terms in small huctuations from the background), (2.16) 


to obtain the held equation at the linearized order: 


= 0 ^ 


SgMN 

_^T^N _ p ^_^2^MN _ ^ 

-VpVqS^^g^^ + V^Sg^^ - 2Rp^q^S^^ 


+RpqS^^g^^ - RS^^ + R^S^^ + R^S^^ + 
+4F^pF^qS^^ + F^S^^ - 2Fl,qS^^g^^] = 0 , ( 2 . 17 ) 


noting that is in the linear order = 0), and we have used the condition Tmn = 0 

for the background held. 
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The generalized Maxwell equation in the linear order in fluctuations takes the form 
1 61' 


— 0 —)■ 2aV mI 


MN 


213V M 2F^pS^^ -2F^ pS 


qPM 


where /mat = Fmn - Fmn- 

The linearized equation obtained from (2.9) is 


-:^MN 


s 


T 


MN 


_ gMN _ S _^ 


1 

2 L 


y2^MAf _ ^ ^ 

-VpVqh^^g^^ + V^hg^^ - 2Rp^Q^h^^ 

FR^^hpqg^^ - Rh^^ + R^h^^ + R^h^^ + Ah 


MN 


+2F^pF^Qh^^ + -F‘^h 


MN 


PpQh'’'^S 


L ‘ 
-MN 


= 0, (2.18) 


-2F^^fp^ - 2F^^fp^ + F^^fpqg^^ = - g^^S. (2.19) 


The important observation is that, using (2.9), the held equation (2.17) becomes 


_^2sMN _ yMyiV^ yMyL^LV ^N^LgLM _ 
(^MiV _ -MNg^ _ 2RpM^NgPq ^ Rp^S^Q-gMN 

FAF^pF^qS^^ + _ 2FlQS^^g^^ = 0 . 


VpVqS^^g^^ + V^Sg^^ 

__l_ 


( 2 . 20 ) 


In the D dimensional Minkowski vacuum, this equation is just the Fierz-Pauli equation 
[25] for the spin-two wave with mass-squared Therefore we recognize that 

corresponds to the massive ghost held. The postulate of no tachyonic ghost requires ^ > 0 
in the present model. 

The linear fluctuation from the background fields can be argued by using above equations. 
In the next section, we consider the solution for spontaneous compactiflcation of spacetime 
as the background. 


III. SPONTANEOUS COMPACTIFICATION IN THE CRITICALLY MODIFIED 
HIGHER ORDER EINSTEIN-MAXWELL THEORY 

In this section, we study the partial compactiflcation of space with a non-trivial flux in our 
model. The extra dimensions are considered to be compactifled into a sufficiently small size. 
A maximally symmetric solution in our model is compactiflcation of the form Mp _2 0 3“^, 
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where Mu _2 is the {D — 2) dimensional Minkowski spacetime and is a two dimensional 
sphere, as in the work of Randjbar-Daemi, Salam and Strathdee (RSS) [26] where D = 6 
is considered. We consider the solution for flat D — 2 dimensional spacetime and study 
its classical stability in this section, though other maximally symmetric spacetimes are also 
interesting. Our notation is that the suffices p, z/,... run over 0,1, 2, D — 3 and the suffices 
m,n,... are used for the extra dimensions D — 2 and D — 1. 

As stated in the previous section, the solution of Tmn = 0 and Jm = 0 is a solution in 
our model. A suitable ansatz for the flux held is 

F = -FMNdx^ A dx^ = -qSmndx'^ A dx^ , (3.1) 

where Smn is the antisymmetric symbol and the constant q is the strength of the two-form 
hux.^ This monopole hux yields the relation 

= \F^9mn ■ (3.2) 

Then, the equation Tmn = 0 can be decomposed to 

~ 2^ F ; (3-3) 

T-mn —F grnn T 9mn T '^^9mn i (3-4) 

for the spacetime of a direct product. As for the hat D — 2 dimensional spacetime, where 
Rfiu = 0, we must take a hne-tuning of parameters as 

F^ = A>0, (3.5) 


and therefore, the solution for the background helds reads 

R = 2F\ R^^ = Ag^^ = 2Fl^, R^^„ = A{5^5^^ - 5^5;) = 2F^^F„ . (3.6) 

Substituting this compactihed background solution into Eq. (2.20), we obtain 


_^2sMN _ y VyL^LM _ 


MN 


^ » (gMiV _ gMKg) ^ g ^ 


(3.7) 


To show complete particle spectrum for the ghost held is rather complicated, but we can 
see that the lowest mass-squared of the ghost is a//? for symmetric tensor helds (which do 


It can be checked that this configuration is the solution of our model. 



not couples to the Ricci curvature of the extra space) and vector bosons of the Kaluza-Klein 
origin (which come from the zero modes of the Lichnerowicz operator on 5'^, — + Rmn)- 

It is also remarkable that the other graviton and vector held huctuations is absent in the 
linearized equation (3.7). 

Thus, taking the analysis by RSS for other huctuation modes [26] into account, in the 
classical and linearized analysis, we hud that the RSS background solution Md -2 is 
stable if q;// 9 > 0, because there is no growing mode.® 

Before closing this section, we note that the effective D — 2 dimensional Newton constant, 
which can be found by the coefficient of the Ricci scalar of the D — 2 dimensional spacetime 
in the action, is not affected by the higher order terms in our model and is independent of 
f3. It might be related to the absence of the graviton and the ghost modes at the linearized 
level mentioned above. 


IV. DE SITTER SOLUTIONS AND COSMOLOGY IN THE MODEL 


In this section, we perform the analysis on the stability of the compactihcation and the 
cosmological solution in our model. We consider an effective potential V for the radius of 
extra space [27] in order to study them. 

We can express the Riemann tensor of the as 

. (4.1) 


where b denotes the radius of the compactihed sphere S^. Then, the equations (3.3,3.4) with 
= 0 impose [26] 


q 


2 


1 



(4.2) 


where the hrst equation is the hne-tuning condition between the parameters. 

The effective potential U as a function of 6, the scale of the extra space, is obtained from 
replacing the background helds by functions of b in —I with 

1 1 

rym rm / Tp2\m ^ rm rm /a o\ 

, [F )n = > ( 4 - 3 ) 


^ Of course, there remain the problem of nonlinear instability and the ghost problem in quantum physics. 
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and i/det gmn oc The effective potential for the RSS model is calculated from —Jq and is 
given as 

(^1 - ^ - a) = -2 + ^ + A^,^ (4.4) 

and then, the effective potential for our model of which Lagrangian is described by (2.7) is 


V{y) = a 


-2 + - + y 

y 


f3A 


AiD - 1 ) 

jy - 1 )^ 

AiD - 1 ) 1/3 


■ , 2D 12iD - 2) 7D- 16 

AiD-2)-Dy + -+ 


y y"^ y^ 


[i7D - 16)f3 + 2iD - A)Py + (4(D - 1)« - Df3)y ‘^], (4.5) 


with y = A6^. A minimum of the potential is found aX y = 1 and attains R(l) = 0. The 
stability of the radius of the extra sphere is determined by the sign of the second derivative 
of the potential: 

V'"(l) = 2 (a + . (4,6) 

If H > 3 and a,/? > 0, the scale b = 1/\/A is found to be stable. We also find that it is 
impossible to create another potential minimum satisfying V = 0 hj any tuning of (A/a in 
the model. 

At a large value of b, V can be indefinitely negative provided that — a > 0. In 

this case, the D dimensional cosmological constant becomes negative. This may not indicate 
instability of compactification against a large fluctuation, because the higher derivative terms 
cause the coupling between the compactification scale and the first and second derivatives 
of an expanding scale factor of the rest of space. 

Therefore, we investigate the cosmological equation including the dynamical scale factor 
and the compactification scale. The metric is usually assumed as 


ds'^ = —dt^ + a^it)dVt\,_^ + b‘^it)dD\ , 


(4.7) 


where bit) denotes the radius of the compactified sphere and dVL\ is the line element an 
the unit sphere. For the maximally symmetric iD — 3) dimensional space, the line element 
is denoted dfl%_^ here, and its Riemann tensor is normalized as R^^ki = kiSldj — Sl6l) with 
k = 0, ±1, where = 1,2,..., D — 3. ait) is the scale factor of — 3 dimensional 

homogeneous and isotropic space. 
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The Riemann curvature is computed with the metric as 



where the dot (') denotes the derivative with respect to the cosmic time t. 

Now we can obtain the effective action of the scale factor a{t) and the compactihcation 
scale b(t) for the RSS model as 




dt a 


D-3i2 


2{D-3)- + 4 + A{D-3)4 + {D-3){D- 4)^^ 


+2 


b^ + 1 
62 


a b 


A64 


-A 


dt a 


D-3i2 


d b 


4{D-3)-- + {D-3){D-4)- 


-a^ + k 


+ 2 - 


a b 

-62 + 1 1 


A 


= Iq , 


(4.9) 


b‘‘ AM 

where the over-all constant which comes from the volume of space has been omitted. The 
action has been integrated by parts after the symbol 

The effective action for our model can be written, by adopting the auxiliary helds, as 


I' 

—fd j dta^~%‘^ 

+2{D-3)s 


2Sr 






a 6 


62 


2A64 


_(fl _ 4)1 - 2 I - 2(fl - 4)11 - + 1 +A 


+4sfe 


.m _ 3)® _ ^ _ m _ 3)®^ _ (-P-3)(T>-4) d^ + k 

a b ab 2 


1 1 
2k¥ ^2^ 


~ 3)s2 + 2sl + [Sq + {D — 3)Sa -|- 2Sfe] I , 


(4.10) 


where we set Sq = sq, Sj = SaS) and S'™ = Sfe(5™. Here, the equations of motion are omitted, 
but they can be obtained from the above action.® 

We seek for the de Sitter (ig'S'^) solution, where k = 0, a(t) = with constant Hq, and 
b = bo =constant, are taken as an ansatz. In this case, all the equations of motion become 


® The Euler equations can be obtained from the action using the command EulerEquations in the add-on 
Calculus' VariationalMethods ‘ for Mathematica@[28]. 
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algebraic equations. With the ansatz, the auxiliary variables are expressed as 

3 - 2A&2 + [D{D - 3)H^ - A]A6^ 

“ 2{D - l)Ab^ ’ 

^ 5-2D + 2{D- 2)Abl - [{D -2){D- 3)Hi + A]Abt 

2{D-l)Abt • ^ ^ 

Substituting these auxiliary relations to the other equations, we obtain the following two 
equations at last: 

-{7D - 16)/3 + 12(A) - 2)/3Abl 

-2 {2{D - l)a + [{D + 2){D- 3){D - 4)H^ + DA]/3} Ab^ 

+4 {2(A) - l)a + [D{D - 3){D - 4)H^ - {D - 2)A]/3} A%1 

+ {(A - 3f{D - 6){D^ -6D + 4:)H*/3 + 2{D - 3){D - 4)H^[2{D - l)a - (A - 2)/3A] 

+ [-4{D-l)a + D/3A]A}A‘^bl = 0, (4.12) 

and 

3(7A - 16)/3 - 24(A - 2)/3Abl 

+2 {2(A - l)a + [(A + 2)(A - 2)(A - 3)H^ + DA]/3} Ab^ 

+ {(A - 3f{D - 2){D^ - 6A + 4)H^/3 + 2(A - 2)(A - 3)H^[2{D - l)a - (A - 2)/3A] 

+ [-4{D-l)a +D/3A]A}A%1 = 0. (4.13) 


The two sets of the solution for b^ and Hq can be found; one is 

(i) = 0 and bl = j, (4.14) 

and another is 

^0 = 3 ^) 3 ^ 

These values for solutions are the same as those of the RSS model [29], as expected (if 
A = 6). In both cases, the values for sq, Sa and Sb are equal to zero. 

Because the simultaneous equations (4.12,4.13) are the fourth order algebraic equations, 
there are other two sets of solutions^ at most. One can hnd, for sufficiently small value for 13, 
only two original solutions (4.14,4.15) are real solutions. For A = 6, only two real solutions 

^ In general, the values for sq, So and Sb are not zero for such solutions. 
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(4.14,4.15) exist for fiKja <~ 2.8. As a model of the universe, the multiple de Sitter phases 
would have some signihcance in cosmology. 

We now study the behavior of the small fluctuation around the above exact solutions. In 
the neighborhood of (i), we set 

6(() = A + AK«). //(«) = - = Ai/((), (4.16) 

VA a 

and Asq = and so on. Then, the linearized equations of motion read 


{D — 3) Asq + 2Asb = 0 , 

(4.17) 

(D — l)Asa + 2(D — 3)A + (D — 1)— As^ = 0 , 

(4.18) 

2{D - l)a [{D - 2)Ab + 2{D - A)AAb 


D - 3 r 


+ [{D -l)iD- 2)a + 4/3A] As, = 0 , 

(4.19) 

{D - l)a [{D - 2)AH - dA^/^Afe 


- [{D -1){D- 2)a - 2{D - 3)/3A] As, = 0 , 

(4.20) 

{D - l)Aso = -2(D - 3)^As, . 

(4.21) 


a 


Note that for /S = 0, the equations give = As;, = Asq = 0, (D — 2)A6 + 2(Zi) — 4)AA6 = 0 
and {D — 2) AH — 4A^/^A6 = 0, which are equivalent to the equations for the cosmology of 
the RSS model [29]. 

The oscillatory behavior of Ab as well as As^ is enough interesting. The variables Asq, 
Asa and Asb oscillate with a different frequency from that of Ab, as seen from (4.18,4.19).® 
Thus, the tuning of initial conditions realizes enough time for accelerating expansion and a 
sufficient inflation can be achieved, though we cannot tell the choice of initial state in the 
present classical analysis of the model. 

Since the analysis of the perturbation from (ii) can be considered and results in hnding 
exponentially growing modes as found in [29], we omit the indication of equations here. 
Some hne tuning gives an enough expansion time for a sufficient inflation as the de Sitter 
phase also in this case. 

® Note that the angular frequency of the oscillation of Aso,a,6 is the square root of in (4.6). 
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V. SPONTANEOUS COMPACTIFICATION IN 
LOVELOCK-HORNDESKI-NONLINEAR MAXWELL THEORY 

In the solution for the spontaneous compactihcation in previous sections, one has found 
the following relation in the background helds: 

Rmnpq = ‘^FmnFpq . (5.1) 

Similarly to the consideration in the previous sections, we come to an idea that the Minkowski 
compactihcation of the RSS model is attained in a higher derivative theory if the relation 
(5.1) is rehected in the higher order correction term. 

We then propose an additional higher order term 

r _ ^sPQRSfpMN opMAfp \(pKL \ 

- ^<Jmnkl[^pq - ^PQ)\^RS~^^ ^RS) 

_ jdMNPQ p a -qMN p I p2 A ( tdMNPQ t? t? a jdMN p2 i p p2\ 

— tt ^KMNPQ — ^-ti KmN + -ti —^yrC ^PmnPRQ — ) 

- HF^ , 

where the generalized Kronecker’s delta is dehned 1 

eiWi 

^N2 

rM2 rM2 
'^Afl (’N2 

r-M„ r-AT® 

^N2 

Now, we write the model 2 by the Lagrangian (therefore, we call the previous model the 
model 1 ) 

L = q.Lq + 7 L 2 • (5-4) 

Note that L 2 consists of the Lanczos-Lovelock Lagrangian (Gauss-Bonnet term) [3-6] and 
the Horndeski’s vector-tensor term [18] and a nonlinear Maxwell term. 

It is easy to see that adding L 2 to the RSS Lagrangian Lq yields no modihcation of 
linearized equation of motion on the RSS background helds after gauge hxing. Therefore, 
the Kaluza-Klein modes remain unchanged from those in the RSS model. It is worth noting 
that the ehective D — 2 dimensional Newton constant in model 2 is not ahected by the value 
of 7 . 

Now, we study the cosmological evolution of the classical background in the model. Al¬ 
though the Kaluza-Klein cosmologies with the Gauss-Bonnet term or the generalized Euler 


rMiM2-"Mp _ 
^NxN2-Np = 


by 


• ■ ■ 

xM2 

■ ■ ■ ^Np 


... 5 


Mp 

Np 


(5.2) 


(5.3) 
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invariant have been studied in many works including [30-35], the case with the RSS-type 
compactihcation with the Gauss-Bonnet term has been studied in a few papers including 
[36],9 

As previously, we characterize the cosmological background as: 

= —dt^ + o?{t)dVt\)_^ + b‘^{t)dVt\ , (5.5) 


and 


F = ^qEmndx'^ A dx” = ^ ^J— Emndx'^ A dx” . 


Then, the Lagrangian L 2 reduces to 


L 2 = 4(11-3)- 
a 


A{D - 4:)-j + {D - 4:){D - 5)- 
a 0 


+ k , h^ + 1 


+ 2 




Ab^ 


b 




+ {D - 3)(Z1 - 4)(dl - 5)(Z1 - 6 ) 


"d^ + k' 


dba^ + k 


a 2 


+8(Zl-3)(Zl-4)(Zl-5)-- 

a b 

a? + k (l)^ + 1 1 




2 +8(Z1-3)(Z1-4)-- 


+A{D - 3){D - A)- 




Ab* 


and we obtain the following action for a and b: 

i^ = J dta^-%‘^\y{D-3){D-A){D-3){D-Q) 
+3{D - 3){D - A){D - 3) 
+A{D-3){D-A) 


1 d^ 2 ^ 

3 aA a? 


1 d^ 6 k db 
3 a? b a? a b 


dH^ + l k-b'^ + l 

+ 




a? 


_4(B-3)(«-4)Z^4-16(Z)-3)A 1 


a b Ab^ } 


(5.6) 


(5.7) 


(5.8) 


a? Ah^ 

Again, we omit exposition of equations of motion here. 

We first explore the de Sitter phase, or the maximally symmetric product spacetime. To 
hud this, we set /c = 0, a{t) = e'^°* and b{t) = do, where ido and 60 sxe constants. We then 


® The compactification with the Horndeski Lagrangian has been considered by one of the present authors 
[ 21 ]. 
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obtain the simultaneous algebraic equations as follows: 

-[a + A{D - 3)(D - A)Hl-i] + 2 [a + 2(D - 3)(D - A)Hl-i\ Abl 
+ [((D -3){D- A)H^ -A)a + {D-3){D-A){D- 5){D - Ab^ = 0 , (5.9) 

a + A{D-2){D-3)H^-f 

+ [((D -2){D- 3)Hl -A)a + {D-2){D-3){D- A){D - Ab^ = 0 (5.10) 

Of course, the solution (i) (4.14) is a trivial solution of above equations, while (ii) (4.15) is 
not a solution for 7 7 ^ 0. One solution of the simultaneous equations (5.9,5.10) with Hq ^ Q 
exists for all positive 7 / 0 ;. The value of 6 o(> :^) increases with the value of 7 /a. The 
existence of the de Sitter phase with large values for Hq and 60 is suitable for the possible 
inflationary phase. 

Next, we draw the phase diagram of the compactihcation scale b. We consider the case 
with k = 0 and D = 6. FIG. 1 shows the phase space spanned by b and b for yA/a = 0,1, 2. 
The small arrows in the hgures indicate the flow in the phase space. One can see that 
6 = 1 /VA is an attractor and the area of attracting region increases with the value of /a}^ 
We conclude that an appropriate choice of the initial condition, which is in the flow passing 
through near by the branch point given by the de Sitter solution of (5.9,5.10), might cause 
a long expansion time. 



FIG. 1. Phase space for model 2. We here take a unit such that A = 1. We show (a) for 7 /a = 0, 
(b) for 7/0 = 1 , (c) for 7/0 = 2 , respectively. 


If ordinary matter is added, the dissipative motion becomes more apparent in an expanding era [37]. 
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VI. SUMMARY AND PROSPECTS 


In this paper, we have considered special cases for extension of Einstein-Maxwell theory 
with higher order terms. The models proposed in the present paper have the same solution 
for spontaneous compactihcation Md_ 2 ^S‘^ utilizing the internal flux as the RSS model. The 
stability of the compactihed vacuum have been conhrmed at the classical level. Nevertheless, 
the models reveal the different cosmological behaviors because of the additional degrees of 
freedom or the existence of the additional de Sitter phases. We should investigate evolutional 
equations by using numerical calculations in future work. 

Our toy models seem to be quite special, but the analysis of the models will exhibit some 
typical behaviors including inflation, and we can regard general higher order corrections as 
a modihcation from the critical models. 

Recently, the black holes in higher derivative gravity attract renewed interest [38, 39]. 
We consider that our model 1 is appropriate to examine magnetized black holes of a novel 
type. 

The subjects we should consider are, for instance: the coupling to matter helds; com- 
pactihcations on [3“^)^ [40] in our models; the inclusion of a p-form flux held such as the 
Freund-Rubin compactihcation [41]; the treatment of our models in quantum cosmology [42- 
44]; the initial singularity problems [45] in the models; the possible supersymmetrization of 
the models; the held content of the model in 3 + 2 dimensions. 

The generalization of our model 1 to more higher order theory is considered straightfor¬ 
wardly. If we introduce an auxiliary held S^, the additional term can be replaced by 


nS 




Ni-N^-iN ‘-’Ml 


qNi qN^—i ( \N PQ'T' 
^Mi' " ^Mr,-i\^M -‘-PQ 


qNl qN„ 

‘^MJ ' ^Ni-Nn ) 


( 6 . 1 ) 


which will be studied extensively elsewhere. 

Finally, we must research the relation to other theories. We suppose that the structure 
of the model 1 owes to the matrix Amnpq, which appears also in the spin-two propagator. 
It would be very interesting if our models were the ehective theory of another theory. 
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